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ABSTRACT 



Through numerical simulations, we study the dissolution timescale of the Ursa Minor cold stellar 
clump, due to the combination of phase-mixing and gravitational encounters with compact dark 
substructures in the halo of Ursa Minor. We compare two scenarios; one where the dark halo is 
made up by a smooth mass distribution of light particles and one where the halo contains 10% of 
its mass in the form of substructures (subhalos). In a smooth halo, the stellar clump survives for a 
Hubble time provided that the dark matter halo has a big core. In contrast, when the point-mass 
dark substructures are added, the clump survives barely for ~ 1.5 Gyr. These results suggest a strong 
test to the A-cold dark matter scenario at dwarf galaxy scale. 

Subject headings: Methods: numerical - Galaxies: dwarf - Galaxies: kinematics and dynamics - 
Cosmology: dark matter 
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1. INTRODUCTION 

The A Cold Dark Matter (ACDM) model has been 
proved to be successful in reproducing structure for- 
mation on large scales. In the standard paradigm, a 
nearly scale invariant power spectrum describes the cos- 
mological primordial density fluctuations. If DM is colli- 
sionless, the clustering of substructure is scale-invariant 
down to the free-strea ming s cale of the CDM parti- 
cle (e.g., EoMameEaDHQOil; [GreeneTaDHoolt). The 
properties of the sub halos in M il ky W a y size galaxies 
have b een studied by iGao et al.l (|2004l ) . iSpringel et all 
(|2008[) and iDiemand et al.l ([20081 ) through collisionless 
DM simulations. They found that 10% of their mass re- 
side in ~ 3 x 10 4 subhalos with masses in the range of 
10 4 — 10 8 M Q , following a CDM subhalo-mass function 
dN/dM oc M" 1 ' 9 . 

At scales of dwarf satellite galaxies, a comparison be- 
tween the model predictions and observations is limited 
by our poor understanding of the baryonic processes in- 
volved in the formation of galaxies. For instance, super- 
nova feedback and gas heating from cosmic sources can 
alleviate the inconsistency between the observed num- 
ber of satellite galaxies in the Milky Way halo and the 
much higher number o f subhalos predicted using A-body 
ACDM simulations dMoore et al.l Il999t Iklypin et af 



2003; 



19991: iBenson et aLll2002HSomervillell2002l : [Ostriker et al 



Font et al.l 120111 ). In this model, many satel- 



lites are either too faint to be detected in surveys or 
dark matter subhalos devoid of a baryonic counte rpart 
(|Bovlan-Kolchin et al.|[2011fc iBovill fc Ricottill2011l ). 

Some observational consequences of the existence of 
substructure within galactic halos have been studied in 
the literature. iRomano-Di'az et al.l ((2008) suggested that 
long-lived stellar bars can be tr i ggere d by a tide from a 
massive subhalo. iKannan et al.l (|2012f ) explored the idea 
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that the interaction of dark matter subhalos with the 
gaseous disks of galaxies may generate density enhance- 
ments in the gas. iJohnston et al.l f|2002T ) studied the dis- 
tribution of carbon stars in the stream of the Sagittarius 
(Sgr) dwarf galaxy and found that the Sgr debris is more 
strongly scattered than would be expected if it would be 
orbiting in a smooth (and thus non-substructured) DM 
halo, but that is entirely consistent with perturbations 
by the Large Magellanic Cloud (LMC) alone. In this 
line of research, iCarlbergl (|2009| ) modeled the interaction 
between streams and halo clumps and concluded that 
stellar streams older than 3 Gyr cannot survive in the 
presence of subhalos with the masses and numbers pre- 
dicted by the ACDM model. Nevertheless, the density 
variations in the star stream detected around the glob- 
ular cluster Palomar 5 and in the galaxy M31 appear 
to be in agreement with the existence of DM subhalos 
(|Yoon et al.ll2011t ICaTlberg et alJl2011[ ). 

The presence of dark-matter subhalos in dwarf 
spheroidal galaxies (dSph) might have observable con- 
sequences and may shed light on the n ature of DM 
partic les. For instance, Pin et all (|2005| ) and ITotanil 
(|2010T l suggested that compact massive dark objects 
in the halo of dSph galaxies could resolve some long- 
standing problem s in the se systems. On the other hand, 
iPenarrubia et all (|2010l) used analytical and A-body 
methods to examine the survival of wide stellar bina- 
ries against encounters with dark subhalos orbiting in 
the DM halos of dwarf galaxies. They found that a large 
fraction of wide binaries can be wiped out due to tidal 
encounters with these dark substructures. The observa- 
tions of large separation binaries would impose a strong 
test to the putative substructure in halos of dwarf galax- 
ies. 

Also lMa et al.1 (j2004T ) studied the gravitational scatter- 
ing effects caused by subhalos on the phase-space distri- 
bution of DM particles in main halos. Their numerical 
experiments indicate that the number and the mass den- 
sity of subhalos could be high enough to flatten the main 
halo's inner cusp within a few dynamical times. 

In this work, we will focus on dSph galaxy Ursa Mi- 
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nor (UMi). UMi has the peculiarity of showing a sec- 
ond stellar density peak (or clum p) which is believe d 
to be a cold long-lived structure (|Palma et al.l I2003D . 
This clump has survived because the underlying grav- 
itational potential in UMi must be close to harmonic, 
which can be accomplishe d if the density pro file of the 
dark halo has a large core dKlevna et al.l(2003t ). In addi- 
tion, iSanchez-Salcedo &: Loral ( 2007t ) established strin- 
gent constraints on the mass and abundance of compact 
objects in UMi in order to preserve the integrity of cold 
small-scale clumps seen in some dSph galaxies. In the 
present work, our aim is to test whether or not the puta- 
tive subhalos with a mass spectrum as that found in colli- 
sionlcss CDM is consistent with the stellar clump found 
in UMi. Our work may shed light on the mechanism 
responsible for the formation of cores in dwarf galaxies 
because some of them predict dissolution of dark-matter 
substructures below a scale of <~ 1 kpc. In fact, the core- 
making mechanisms proposed so far fall in three broad 
categories: supernova feedback (e.g., Mo & Mao 2004; 
Mashchenko et al. 2006; Governato et al. 2010, 2012; 
Pontzen & Governato 2012; Maccio ct al. 2012), dynam- 
ical friction from infalling baryonic clumps (El-Zant et 
al. 2001; Romano-Diaz et al. 2009; Goerdt et al. 2010), 
or a change in any of the basic properties of dark mat- 
ter particles [e.g., collisional, annihilating, decaying or 
warm dark matter] (Spcrgcl & Stcinhardt 2000; Kapling- 
hat, Knox & Turner 2000; Cen 2000; Sanchez-Salcedo 
2003; Avila- Reese et al. 2001). Models in the latter cat- 
egory are expected to help erase dark-matter subhalos 
in dwarf galaxies, e.g., by mass stripping if dark mat- 
ter is collisional or by suppression of the power in the 
mass spectrum at small scales if dark matter is warm, 
whereas models in the first two categories are expected 
to preserve them. 

This article is organized as follows. In fJ21 wc describe 
some properties of UMi and its clump. The initial con- 
ditions for our A-body simulations are described in 
The results of the simulations are shown in SJD Finally, 
we present our conclusions in SJ3J 

2. UMI AND ITS CLUMP 

UMi is a dSph galaxy satellite of the Milky Way, lo- 
cated at a Ga l actoc entric distance of Rgc = 69 ± 4 kpc 
(|Grebel et al.l I2003D . Dynamical studies suggest that 
UMi has a mass-to-l i ght r atio larger than 60Af Q /L Q 
(e.g. IWilkinson et al.l 120041 ) . Rescaling the results of 



IGilmore et al.l (12007ft and taking a luminosity of L\ 



1.1 x 10 6 L^ (jPalma et al.ll2003D . the total mass within 
0.6 kpc is 6.3 x 1O 7 M , which results in a mass-to-light 
ratio M/L v > 60M Q /L Q for UMi. It has to be noted 
that if we take L v = 3 x 10 5 i Q (IGrebel et al.l I2003D 
then we obtain a M/Ly ~ 2OOM Q /L0 (jGilmore et al.l 
|2007|) . Deeper data for dwarf galaxies often reveal a 
larger angular extent of their stellar compon ent (e. g., 
lOdenkirchen et al.l 120011: iKniazev etail I2009D ; neverthe- 
less, we adopt the lower luminosity value for UMi. The 
high value of M/L in UMi implies that it is one of the 
most dark matter dominated dSph galaxies in the Local 
Group. 

UMi's stellar King core radi us along the semima jor axis 
is 17.9 arcmin (~ 0.4 kpc) (jPalma et al.ll2003lh UMi 
shows a large ellipticity in the shape of the inner isoden- 



sity contours of th e surface density of stars (e ~ 0.54) 
(jPalma et al.l 120031 ). but the most remarkable feature in 
UMi's structure is t he second off-centered density peak 
(jKlevna et al.|[l998l ). One of these two peaks is located 
on the north-eastern side of the major axis of UMi at 
a distance of ~ 0.4 kpc from UMi's center. The radial 
velocity distribution of the stars in UMi is well fitted 
by two Gaussians, one representing the underlying back- 
ground (8.8 km s _1 ) and the other representing the ve- 
locity dispersion of the second peak (0.5 km s _1 ). The 
stars in the vicinity of this peak comprise a kinematically 
distinct cold subpopulation: a dynamically cold stel- 
lar clump. The most appealing interp retation of UMi's 
clump is that it is a disrupted cluster (jRead et al.ll2006l) 
that orbits in the plane of the sky, and that has survived 
in phase-space becau se the underlying g ravitational po- 
tential is harmonic ([Klevna et al.l 12003). This implies 
that the dark halo density profile in UMi should have a 
core (and not a cuspy profile as predicted by the ACDM 
par adigm) and t hat this core must b e larg e (~ 500 pc, 
see IKlevna et al]|2003l and ILora et al.ll2009f) . 

The stellar population of the UMi dwarf galaxy is very 
old with an age of 10-13 Gyr; virtually all the stars were 
form ed 10 Gyr ago, and 90% of them were formed 13 Gyr 
ago (jCarrera et al.l I2002t IRocha et all 1201 ID . This im- 
plies that star formation was halted at a redshift of z ~ 2, 
probably because supernova explosions were able to re- 
move all the gas. In fact, deep observa tions have not 
detected gas in UMi (jYound [19991 l2000lh Gas replen- 
ishment from dying stars and stellar winds is expected 
to be very small. If the gas returned to the interstellar 
medium by intermediate- and low-mass dying stars were 
distributed across the galaxy as the stars are, the central 
gas density wo uld reach a value of ~ 4.4 x 10~ 3 M Q pc~ 3 
after 10 Gyr (C arrcra et al.l I2002D . This corresponds 
to a particle density of 5 x 10~ 3 cm~ 3 , much smaller 
that the typical lower threshold s for star-forming regions 
(~ 1000 cm" 3 ) (|Shu et al.lll987t) . Since we are interested 
here in the dynamical evolution of the stellar clump over 
the last 10 Gyr, once essentially all the gas has been 
removed, we can discard both gas processes and stellar 
evolution in our simulations. 

3. JV-BODY: INITIAL CONDITIONS AND THE CODE 

3.1. The dark matter component of UMi 

We performed A-body simulations of the UMi dSph 
galaxy, which include a baryonic component (the bulk 
stellar component and the stellar clump) embedded in a 
spherical (live) dark matter halo. The dark matter halo 
density profile selected is defined as 



p{r) 



Pi) 



(1) 



Here po accounts for the central density in the case 
of 7 = and r s is the scale length. A NFW profile 
(|Navarro et al.l fl996() is obtained for (a, /?, 7)=(1, 3, 1). 
We have chosen the values (a, fl, 7)=(1, 3, 0) in order to 
have a cored DM mass density profile. 

It has been found that cored profiles (and not 
cuspy profil es) are in b e tter ag reement wi t h UM i 
dynamics (jKlevna et al.l 120031: ILora et al.l I2009D . 
Some other dSph galaxies in the Local Group such 
as Fornax, Sculptor, Carina, Leo I and Leo II arc 
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also believed to p ossess cored dark halos instead 
: cu spy profiles felevna et al.l^j2003j ; iGoerdt et al 



2006L [2010t ISanchez-Salcedo eta. 



iGilmore et al 



20071 iBattadia et all 120081: lAmorisco et all 1 2011 



Walker fc Penarrubial |201H fjardel fc Gebhardtl 12012 : 
Agnello fe Evansll2012| ). Moreover, it is found from high 
resolution observations of the rotation curves for DM 
dominated low surface brightness (LSB) galaxies that 
DM halos h ave density profiles compatible with flat 
central cores (|de Blok et al.l[200ll : IChen et al.lf2(TToh . 

We explored two different cored halos; a small-core 
halo with a scale length of 0.91 kpc (which corresponds 
to a core radius of ~ 0.4 kpc) with a total mass of 
M = 2 X 1O 9 M , and a big-core halo with a scale 
length of 2.2 kpc (core radius ~ 1 kpc) and a total 
mass of M = 3 x 10 10 M Q . For our halo models we 
have a mass within a radius of 0.39 kpc (the clump's 
orbit) of ~ 10 7 A/ Q for the small-core halo, and a mass 
of ~ 3 x 10 7 M Q for the bi g-core halo mo del. This is 
consistent with the results of Strigai Tet al.l (2007, 2008) 
who considered dark halos compatible with the observed 
stellar kinematics of the classical dSph galaxies, includ- 
ing the dwarf galaxy UMi. They found that, for realistic 
density profiles, the mass interior to 300 pc is ~ 1O 7 M 
for all dSph galaxies in the Milky Way halo. 

To generate the initial conditions of the dark matter 
particles of the smooth background dist ribution, we use d 
the distribution function proposed by iWidrowl (|2000ft . 
and the velocity dispersion of the system was taken to 
be isotropic. 

In order to explore the dynamical effects of a clumpy 
halo, we generated the mass spectrum of substructures 
in the halo following the power-law mass distribution 



dN 
dM 



M 
Mo 



(2) 



(IGao et all 120041: ISpringel et ail 120081: IDiemand et~aH 
12008ft . The values of ao and M are fixed once the lower 
M m i n and upper M max mass limits of the substructure 
particles and the fraction (/) of the halo mass comprised 
in substructures, are given. 

The mass of dark matter particles used in our simula- 
tions without subhalos is 2 x 10 3 M Q , which corresponds 
to our maximal mass resolution for the halo particles. 
Then it is a natural choice to take M m i„ = 10 3,3 Mq 
as the lower mass limit. Adopting the subhalo mass 
function found in CDM simulations, the subhalos are ex- 
pected to have masses up to M max < 0.1M wr , where 
M vir is the virial mass of the halo (IPenarrubia et all 
2010!). Given that the estimated average virial mass 



of the dSph galaxies is of the order of M„ 



10 9 M,r 



(jPenarrubia et al.l 120081: IWalker et all 12009ft . the upper 
mass limit could be as large as 10 8 M Q . It then be- 
comes reasonable to take the generous upper mass limit 
as M max ~ 10 7 M p . 

Equation ^ with /, M m i n and M max known, deter- 
mines the number N$ub of substructure particles. In 
order to construct a clumpy DM halo, we choose from 
the original smooth DM halo, a number Nsub of ran- 
dom particles and replace their masses according to the 
power- law in Equation 



3.2. The baryonic component of UMi 
The spatial stellar density profiles of elliptical systems 



are commonly described as a power law in radius 



(3 - 7 )M* 



-in 



(3) 



where M* is the total stellar mass and a is the scale 
radius. For f3 = 4, these models have s imple analytic 
properties and are ca lled Dehnen models (lDehnenlll993l : 
iTremaine et al.l [1994ft . In these models the densities are 
proportional to r -4 at large radii and diverge in the cen- 
ter as r -7 . To model the bulk stellar co mponent in UMi 
we used a Dehnen model (|Dehnenlll993[) where 7 = 3/2. 
This slope most closely resembles the de Vaucouleurs 
model in surface density. We took the scale radius for 
the underlying s tellar component in the UMi galaxy t o 
be a = 0.4 kpc (jKlevna et al.l 119981: iPalma et all 120031) . 
We set the total mass of the stellar component to be 
M* = 9 x 10 5 M Q taking the typical value of the mass- 
to- light ratio T* = 3. 

We performed an iV-body simulation with the DM halo 
and the underlying stellar components together. The 
resulting system was found to be stationary for a Hubble 
time (i.e. the density profile of the stellar component and 
the velocity dispersion stayed approximately constant for 
a Hubble time). The stellar velocity dispersion (er ~ 
10 km s _1 ) was in agreement with the velocity dispersion 
of the stellar component found in UMi (a = 9.3, 12, 9.5 ± 
1.2 km s" 1 : iWilkinson et al.l [200l [Gilmore et all 120071: 
IWalker et al.1 120091) . 

Finally, for the initial density profile of the stellar 
clump, we take 



Pc(r) = p exp(-r 2 /2r^), 



(4) 



with the clump ra dius r c between 12 pc and 35 pc 
(jParma et al.l 12003ft . The clump's velocity dispersion was 
set to 1 km s _1 . The clump was dropped at a galactocen- 
tric distance o f 0.39 kpc in a circular orbit in the (x, y) 
plane (see also San chez-Salcedo fc Lora II2010D . Since the 
clump's stars have the s ame color as the un derlying stel- 
lar population in UMi ([Klevna et al.l 11998ft . we assume 
that the U-band mass-to-light ratio M /Ly of the clump 
is the same as it is for the underlying stellar compo- 
nent (say, T* = 3). Thus, the mass of the clump M c is 
~4 x 10 4 M o . 

3.3. The code 

Since the internal two-body relaxation timescales for 
the three components (clump, underlying stellar com- 
ponent and halo) arc much larger than one Hubble 
time, this system can b e represented as collisionless 
(jBinnev fc Tremaind 12008ft . We simulated the evolution 
of the UMi dwarf galaxy (stellar clump, underlying stel- 
lar compo nent and DM halo) u sing the A^-body code su- 
perbox (|Fellhauer et al.l [20001 ) . superbox is a highly 
efficient particle-mesh, collisionless-dynamics code with 
high resolution sub-grids. The parameters of each model 
are given in Table [T] 

In our case, superbox uses three nested grids centered 
in the center of density of the UMi dSph galaxy. We used 
128 3 cubic cells for each of the grids. The inner grid is 
meant to resolve the inner region of UMi and the outer 
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grid (with a radii of 1000 kpc for all cases) resolves the 
stars that are stripped away from UMi's potential. The 
spatial resolution is determined by the number of grid 
cells per dimension (N c ) and the grid radius (r gr id)- Then 
the side length of one grid cell is defined as I = • To 
study the convergence of the results with different res- 
olutions, we performed three simulations: a simulation 
with 128 3 cubic cells, one with higher (256 3 cubic cells) 
and one with lower (64 3 cubic cells) resolution of model 
M3 (see model parameters in Table [TJ. We obtained 
practically the same results for the three resolutions. We 
conclude that 128 3 cubic cells suffices to achieve a robust 
accuracy. For N c = 128, the resolution, which is of the 
order of the typical distance between the particles in the 
simulation, is given at the two last columns of Table O 

SUPERBOX integrates the equations of motion with a 
leap-frog algorithm, and a constant time step dt. We se- 
lected a time step of = 0.1 Myr in our simulations in 
order to guarantee that the energy (for the isolated com- 
ponents) is conserved better than 1%. The properties of 
the three modeled components of UMi (DM halo, under- 
lying stellar component and stellar clump) used in the 
simulation are shown in Table [2j Also, the radii of the 
inner, middle and outer grids for each of the components 
are given. 

4. RESULTS 

Our iV-body simulations were carried out from a time 
t = to t = 10 Gyr using the code described in Sec- 
tion 13.31 In Table Q] we list the parameters used in the 
simulations (see also tj3.2[) . For the models Ml — A/4 
(small-core halo) we used 10 6 particles to represent the 
halo component, 10 5 particles for the bulk stellar com- 
ponent, and 10 4 particles for the stellar clump. For the 
models M5 — A/8 (big core halo) we used the same par- 
ticle number mentioned before for the underlying stel- 
lar component and for the stellar clump, but we used 
1.5 x 10 7 particles for the halo component. This was 
done in order to have the same mass per particle in both 
small (2 x 10 9 M o ) and big (3 x 10 10 M Q ) core DM halos. 
In total we have a set of 8 simulations (see M1-M8 in 
Table [TJ. 

4.1. Halo without substructure 

In Figure [T] we show the evolution of the clump in 
the models Ml (clump radius r c = 12 pc) and M3 
(clump radius r c = 35 pc) for a halo with a small core 
(r s = 0.91 kpc). At t = 1 Gyr in model Ml, the clump 
inflated its radius by a factor ~ 1.5 due to gravitational 
encounters with the halo particles, but it preserves its 
identity (see panel (b) of Figure [T]). In model M3, the 
clump appears to be very disrupted already at t = 1 Gyr 
(see panel (e) in Figure []}. We can say that, regard- 
less the clump radius, the clump is totally destroyed at 
t = 2.5 Gyr in both models (see panels (c) and (/) in 
Figure [J) . 

In order to quantify the destruction time of UMi's stel- 
lar clump in our simulations, we built a map of the sur- 
face density of the stellar clump in the (x, y)-plane at 
any given time t in the simulation. We sample this two- 
dimensional map searching for the 10 x 10 pc size parcel 
that contains the highest mass (number of clump par- 
ticles). This region will be centered at the remnant of 



the clump. We define that a clump is destroyed when 
this region has reached a density of 1M pc -2 . When 
such a small value of the surface mass density is reached 
by UMi's clump, the column density of the clump is so 
low that it would be unrecognisable from the underlying 
stellar component, and would be thus undetectable. 

Figure [5] shows the surface density of mass as a func- 
tion of time, where the black line represents the surface 
density at which the clump is disrupted. The destruc- 
tion times for both clumps in the small-core case are 1.8 
and 0.9 Gyr for r c = 12 pc and r c = 35 pc, respectively 
(see last column of Table [T]) . This implies that a halo 
core of ~ 0.4 kpc is not large enough to guarantee the 
survival of th e clump . This result is in agreement with 
iKlevna et all 's (|2003f) statement that if the core radius 
of the DM halo is equal to or smaller than the clump's 
orbit, the clump will get destroyed within ~ 2 Gyr. We 
recall that, in this case, the DM halo has a core radius 
of ~ 0.4 kpc and the clump has a galactocentric distance 
of 0.39 kpc. 

In Figure [3] we show the evolution of the clump for 
both r c = 12 pc (A/5) and r c = 35 pc (A/7), embedded 
in a halo with a large core (r s = 2.2 kpc). In both 
cases the clump survives for ~ Hubble time. The small 
(r c = 12 pc) clump expands up to ~ 2 times of its original 
size in the first 5 Gyr (sec panel (b) of Figure [3]) . After 
that, it maintains this size over one Hubble time. In the 
M7 model (r c = 35 pc), the clump loses some particles 
(sec panel (e) of Figure [3]) and slightly reduces its initial 
size (see panel (e) of Figure^. Then the clump continues 
losing some particles and, as a consequence, it shrinks 
its initial size by a factor of ~ 2.5 by the end of the 
simulation (~ Hubble time), but remains undestroyed. 
The survival of the clump can also be seen in Figure [2j 
the surface mass densities in the M5 and Ml models lie 
well above the clump's destruction line (indicated by the 
horizontal black line). 

This set of simulations tells us that a halo with a large 
core allows the survival of both clumps for ~ Hubble 
time. This result reaffirms the belief that t he UMi dSph 
galax y should have a large core DM halo (jKlevna et al.1 
2003;), instead of a cuspy DM density profile. 

4.2. Halo with substructure 

In this section we present the same simulations as those 
in H4.ll but including substructure in the DM halo with 
/ = 0.1, M max = 10 7 A/ Q and M min = 1O 3 3 M (which is 
the maximum mass resolution for the DM halo particles). 
For the halo with a small core (where the total mass of 
UMi is A/ = 2 x 10 9 A/ Q ), the number of substructure 
particles is of N su t, = 5836. For the halo with a large 
core (M = 3x 10 10 M Q ), we require N sub = 22483. 

In Figure [4] we show snapshots of the models M2 and 
A/4 (small-core clumpy DM halo) at t — 0, 1 and 2.5 Gyr. 
The purple circles represent the DM particles, the small 
light gray circles represent UMi's extended stellar com- 
ponent particles, the black points represent the particles 
that make up the stellar clump and in green we show the 
substructure particles. From Figure [4] (see panels a, b 
and c), we can see that the small (r c = 12 pc) clump has 
increased its size by a factor of ~ 2 after 1 Gyr, slightly 
larger compared to the case without substructure. We 
built a map of the surface density of the stars initially 
in the clump (see Figure [5]), similar as that in the case 
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without substructure. A clump with an initial radius 
of r c = 35 pc survives ~ 1 Gyr, which is also approxi- 
mately the same time as in the non-substructured case 
(see last column of Table Q} . These results imply that 
the dissolution by tidal forces in a DM halo with a small 
core plays a major role in the dynamical evolution of the 
clump, more than the interaction between the massive 
substructure particles and the particles of the clump. 

In Figure [6] we show the number of substructure par- 
ticles with masses > 1O 4 M inside a sphere of 0.39-kpc 
radius for models M2 and MA. This radius corresponds 
to the radius at which the clump is orbiting, and thus 
massive subhalos within that radius will have a major 
effect in the clump's destruction. In these models, the 
number of particles more massive than 10 4 Mq, inside 
the clump's orbit, ranges from 1 to 9, with a mean num- 
ber of 5 particles. 

In Figure [7] wc show snapshots for the clumpy halo with 
a big core (models M6 and MS) at t = 0, 0.5 and 1.5 Gyr. 
In both simulations, the clump appears to be enlarged 
at t = 0.5 Gyr (see panels b and e of Figure [7]) and, at 
t = 1.5 Gyr, it is practically disrupted. The clump is 
dissolved at 1.6 Gyr in model M6, and at 1.4 Gyr in 
model M8. The number of substructure particles with 
mass larger than 10 4 M Q is shown in Figurc[5] In the case 
of a halo with a big core, the number of particles with 
masses > 10 4 Af Q inside a 0.39 kpc sphere ranges between 
8 and 23, with a mean number of 16. In the models with 
a large core, but without substructure, the clumps (12 
and 35 pc) remain undestroyed for ~ Hubble time, which 
means that the destruction of the clump in the clumpy 
halo case is due to the random walk in momentum space 
that the stars in the clump undergo by the collisions with 
the substructure particles. This effect is so important 
that even with a core as large as ~ 1 kpc, the clump 
does not manage to survive the continuous encounters 
with the massive substructure particles of the DM halo. 

4.3. A comparison with analytical estimates 

iSanchez-Salcedo fc Lora I (|2007l ) studied the abun- 
dance of Very Massive Objects (VMO) in the DM halo 
of Fornax and UMi. These VMOs can be compared with 
the subhalos that we have studied in this paper for UMi. 
In the impulse approximation, they found that, if the 
progenitor cluster became unbound immediately after 
formation, the mass of the VMOs should be 

M h < 2.5 x 1O 3 M ( fPdm 1 ( gdm , ) 

° V0.1M Q pc- 3 y V20kms-V 

(5) 

in order to have a clump as dynamically cold as observed. 
Here pdm and o~dm are the density and velocity disper- 
sion of the dark matter particles in the UMi halo. For 
/ = 0.1, this upper limit implies Mh < 2.5 x 10 4 M Q . 
Therefore, the maximum number of VMO within the or- 
bital radius r or b is ~ 47r/r , rb p dm /(3M/ l ) ~ 100. In the 
case of a mass spectrum, the maximum number of sub- 
halos with masses < 2.5 x 10 4 M Q are expected to be 
smaller because the dissolution time of the clump due to 
collisions with subhalos of mass M goes as oc M 2 dN(M). 



Our TV-body simulations have shown that, if the subhalos 
follow a mass spectrum and if self-gravity of the clump 
is taken into account, ~ 16 subhalos with masses larger 
than 10 4 M Q are enough to disrupt the clump if the tidal 
perturbation by the dark halo is also included. 

5. CONCLUDING REMARKS 

We studied the dynamical consequences of the putative 
substructure in UMi's DM halo in the stellar component. 
We ran iV-body simulations of a stellar clump orbiting 
in a live cored DM halo. If the distribution of mass in 
the dark halo is smooth, a dark halo with a scale length 
of 0.91 kpc, which corresponds to a core of ~ 0.4 kpc, 
cannot preserve the integrity of the clump. On the other 
hand, for a dark halo with a big core (scale length of 
2.2 kpc, core radius of ~ 1 kp c), the clump survi ves for 
approximately a Hubble time ([Klevna et al.|[2003T ). 

When 10% of the original DM halo total mass resides 
in compact subhalos, the clump dissolves in roughly the 
same timescale in the small core case. This means that 
the dissolution effects by tidal forces by the small-core 
DM halo are greater than due to the gravitational scat- 
tering with the massive substructure particles. 

In the big-core case there is a large number of particles 
in substructure (N su b = 22483) with a mean number 
of 16 particles inside a sphere of radius 0.39 kpc (the 
radius of the orbit of the clump) with masses greater 
than 10 4 M Q throughout the simulation. The effect of 
the substructures over the clump in this case results in 
the complete destruction of the clump within ~ 1.6 Gyr 
for the clump with r c = 12 pc and within ~* 1.4 Gyr for 
the clump with r c = 35 pc. In a smooth dark-matter 
halo, the large DM core halo ensures the longevity of 
the clump for almost one Hubble time, but in a clumpy 
halo, the clump is erased because of the random walk in 
momentum space that the stars in the clump undergo by 
collisions with the massive substructure particles. 

It remains to study a more realistic scenario where the 
subhalos are not point particles (in which case resembles 
more the case of VMOs) but are extended perturbers. 
It would be worthwhile to quantify the effect that these 
extended subhalos would imprint compared to the point 
mass case studied here, and see in a more general case 
if subhalos in dSph galaxies impose a strong test to the 
ACDM scenario. Preliminary simulations suggest that 
the disruption timescale of the clump increases by 60%- 
70% if subhalos are extended, which is not enough to 
account for the survival of the clump. A thorough study 
will be presented elsewhere. 
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Fig. 1. — Evolution of a clump embedded in a smooth halo with a small core. The time is given at the bottom right corner. The top 
panels (a, b and c) show the case where the radius of the clump is r c = 12 pc (model Ml). The bottom panels (d, e and /) show the case 
where the radius of the clump is r c = 35 pc (model M3). The purple circles represent the soft halo DM particles; the small light gray 
circles indicate the extended stellar component; and the black dots are for the particles of the stellar clump. 




Fig. 2. — Surface density of the mass of UMi's clump map in the (x, i/)-plane at any given time t in the simulation for the four models 
without substructure (see the models Ml, M3, M5 and Ml in Table[l}. The black line shows the destruction (IMq pc - 2 ) line. 
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Fig. 4. — Snapshots at t = 0, 1 and 2.5 Gyr for the small-core DM halo with 10% of the DM halo's total mass in substructure (green 
circles). The top panels show model Ml and the bottom panels show model MA. The color code of the particles is the same as in FigurefTI 
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Fig. 5. — Same as Figurc[2]for the four models with dark-matter subhalos (see models M2, MA, M6 and M8 in Table[T} 
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Fig. 6. — Number of halo particles with mass > 10 4 Mq, inside a 0.39 kpc sphere around the center of density of the simulation (i 
inside the clump's orbit) for model M2 (black line) for model MA (red line). The number of particles was computed at each 0.1 Gyr. 
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Fig. 7. — Same as Figurc[4]but at times t = 0, 0.5 and 1.5 Gyr. The top panels show model M6 and the bottom panels show model M8. 
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TABLE 1 

Parameters used in the simulations. Here r a accounts for the DM halo scale length; M is the total mass of the halo, and 

r c is the radius of the stellar clump in UMi. 



Model 


sub- 




M 


r c 


destruction time 




structure 


[kpc] 


[M@] 


[pc] 


[Gyr] 


Ml 


no 


0.91 


2 x 10 9 


12 


~ 1.8 


M2 


yes 


0.91 


2 x 10 9 


12 


~ 2.1 


M3 


no 


0.91 


2 x 10 9 


35 


~ 0.9 


Mi 


yes 


0.91 


2 x 10 9 


35 


~ 1 


M5 


no 


2.2 


3 x 10 10 


12 


> 10 


M6 


yes 


2.2 


3 x 10 10 


12 


~ 1.6 


M7 


no 


2.2 


3 x 10 10 


35 


> 10 


M8 


yes 


2.2 


3 x 10 10 


35 


~ 1.4 
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TABLE 2 

Setup properties of the particles in the simulations. Total mass, mass per particle and number of particles of the three 
UMi modeled components, the radii of the grids and the resolution I of the middle grid (for N c = 128) ARE SHOWN. 



Component 


M 


m p 


N 


inner grid 


middle grid 


'l28 




[Mq] 


[Mq] 


10 6 


[kpc] 


[kpc] 


[kpc] 


Halo with small core 


2 x 10 y 


2 x 10 J 


1 


1 


20 


0.322 


Halo with large core 


3 x 10 10 


2 x 10 3 


15 


2 


30 


0.483 


Extended stellar component 


9 x 10 5 


9 


0.1 


1 


10 


0.161 


Clump with r c = 12 pc 


4.03 x 10 4 


4.03 


0.01 


0.012 


0.12 


0.0019 


Clump with r c = 35 pc 


4.03 x 10 4 


4.03 


0.01 


0.035 


0.35 


0.0056 



